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Abstract 

We calculate within a Bose-Hubbard tight-binding model the matter-wave flow 
driven by a constant force through a Bose-Einstein condensate of 87 Rb atoms in var- 
ious types of quasi-onedimensional arrays of potential wells. Interference patterns 
are obtained when beam splitting is induced by creating energy minigaps either 
through period doubling or through quasi-periodicity governed by the Fibonacci 
series. The generation of such condensate modulations by means of optical-laser 
structures is also discussed. 
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1 Introduction 



A Bose-Einstein condensate (BEC) is a gas in which a macroscopic number 
of massive particles reside in the same quantum state (see [1] for a review of 
work done on quasi-pure BEC's produced since 1995). Experiments aimed at 
revealing the coherence of a BEC have demonstrated its matter-wave prop- 
erties. In particular, condensate interferometry can be realized by splitting a 
BEC into two parts with a definite phase relationship, these parts being then 
brought into overlap and interference as for an optical laser beam that has 
gone through a beam splitter. Coherent splitting of a BEC has been achieved 
by optically induced Bragg diffraction [2] and a number of ingenious methods 
have been devised to extract a collimated beam of atoms from a BEC (see 
e.g. [3]). 

A quasi-onedimensional (ID) array of potential wells is created for an atomic 
BEC by the interference of two optical laser beams which counterpropagate 
along the z axis, say, and are superposed on a highly elongated magnetic trap. 
Such an optical lattice provides an almost ideal periodic potential and has 
allowed the study of Bloch and Josephson-like oscillations [4,5,6] and of the 
mechanisms by which decoherence arises as in the transition from a superfluid 
to a Mott-insulator state [7]. 

The ID optical lattice can be modified by means of auxiliary laser beams [8]. 
In particular, its periodicity can be doubled by adding two beams that are 
rotated by angles of 60° and 120° with respect to the z axis. For a suitable 
choice of the phases the potential seen by the BEC atoms takes the shape 

U(z) = U [sin 2 (kz) + (3 2 sin 2 (kz/2)} , (1) 

where U is the potential well depth, f3 2 is the relative energy difference be- 
tween adjacent wells, and k is the laser wavenumber determining the distance 
d of adjacent wells as d = ir/k. In solid-state terminology, the doubling of the 
period when f3 2 ^ causes the opening of a minigap in the energy spectrum as 
a function of f3 2 . A BEC driven through such a lattice by a constant force is co- 
herently split by a combination of Bragg diffraction and of tunnelling through 
the minigap. In steady state a prominent interference pattern is produced and 
can be observed by monitoring the outgoing particle flow [8]. 

In this work we show how a quasi-periodic Fibonacci array of potential wells 
could be created by optical means and evaluate its density-of-states structure 
to display a series of approximate minigaps. We then show that this structure, 
unlike a simple periodic structure but similarly to a period-doubled structure, 
leads to an interference pattern under steady-state drive of the BEC by a 
constant force. The model and the behaviours of periodic arrays are briefly 
reviewed in Sec. 2 and Sec. 3, respectively. The Fibonacci array is treated in 
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Sec. 4, while Sec. 5 offers some concluding remarks. 



2 The model 



We use a ID tight-binding Hamiltonian for the BEC atoms and a Green's 
function approach to evaluate their linear transport coefficient through the 
array of potential wells [8]. The Bose-Hubbard Hamiltonian for N bosons 
distributed inside n s wells is 

^ = X;[^|i><i|+7i(IO<< + l| + l< + l><<l)] • (2) 

Here, the parameters Ei and 7, represent site energies and hopping energies, 
respectively, and depend on the number of bosons in the well labelled by the 
index %. In a tight-binding scheme the ID condensate wavefunction in the 
i-th well is a Wannier function for the bosons in the external potential and, 
according to the early work of Slater [9], can be written as a Gaussian function 
of longitudinal width a z , 

<f>i(z) = h( Zi ) exp[-(z - z i ) 2 /(2a 2 z )} . (3) 

The parameters entering the effective Hamiltonian are given by 

h 2 v 2 1 



Ei = J dz<pi(z) 



2m v ' 2 
for the site energies and by 

h 2 v 2 , , 1 



+ U(z) + -g bb \<P % (z)\ 2 -maz + C 



Hz) (4) 



li = J dz<f)i(z) 
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2m w 2 



<t>i+i{z)- (5) 



for the hopping energies. In Eqs. (4) and (5) m is the boson mass, a = F/m 
is the acceleration due to a constant external force F acting on the bosons, 
g bb is the effective ID boson-boson interaction parameter, and C is a constant 
accounting for transverse effects in a cigar-shaped trap (for the determination 
of the parameters see Ref. [8]). We remark that in a tight-binding approach 
nonlinear interaction effects enter the self-consistent calculation of the axial 
width a z , resulting in a broadening of the Gaussian function, and also modify 
the on-site energies and the hopping energies. This approach is justified in the 
case of weak boson-boson coupling as for a 87 Rb BEC, on which we focus in 
this paper, and should be improved in a strong-coupling situation as is met 
on the approach to a Feshbach resonance [10]. 

In the Green's function method the calculation of the transmittivity of bosonic 
matter waves through the array of potential wells does not require an ex- 
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plicit solution of the Hamiltonian (2). The array is reduced by a renormaliza- 
tion/decimation technique to a single "dimer" [11], to which an incoming lead 
and an outgoing lead are connected. The steady-state transport coefficient is 
inferred from the scattered wavefunction of the leads in the presence of the 
dimer. 



Period-doubling of the array as described by the expression of U(z) in Eq. 
(1) yields an interference pattern as a function of the ratio T B /r, where 
T B = hk/ma is the period of Bloch oscillations and r is the average time 
needed for tunnelling twice across the minigap. The array acts in this case as 
an interferometer for bosonic matter waves and we construct below its optical 
analog. In essence the minigap plays the role of a medium with large refrac- 
tive index, across which an evanescent wave couples two layers that allow real 
wave propagation. In the case of a Fibonacci array, on the other hand, peri- 
odicity is lost but quasi-periodicity induces the opening of a number of rather 
sharp depressions in the density of states ( "quasi- minigaps"). The resulting 
fragmentation of the spectral density modifies the interference pattern, but 
does not erase it. 



3 Interference from period doubling 



When f3 2 = the ideal infinite lattice generated by U(z) has a single period 
and its low-energy spectrum is that of a one-band system, as is shown in the 
left panel of Fig. 1. Period doubling causes the opening of a minigap in the 
total density of states (DOS) (see right panel in Fig. 1). The calculation of the 
DOS has been carried out by recursive algorithms as those in Refs. [13,14] and 
full details will be given in a later publication. The energy width AE of the 
minigap is fixed by the energy difference \Ei — E i+1 \ between two adjacent sites 
and in the limit of an infinite lattice no states are present inside the minigap. 
In real systems as in the experiments at LENS (see for example [12]) the 
condensate occupies about 100-200 wells. As a result of finite-size effects the 
gap is not completely empty and the bosons can easily be transferred by 
tunnelling between the two sub-bands. 
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Fig. 1. Total DOS of a very long lattice (1000 sites) with a single period (left panel) 
and a doubled period (right panel), as a function of energy E referred to the central 
energy Eq and with 4i being the total spectral width. In the single-period lattice Eq 
is the site energy Ei and t the hopping energy ji, while in the double-period lattice 
E = (Ei + E i+1 )/2 and t = [(Ei - E t+1 ) 2 /2 + 4 7 ?] 1 /a/2. 



The addition of a potential maz causes a tilt of the bands in space and the 
density of states depends on both position and energy. Thus the condensed 
bosons are driven through the lattice and explore the whole band (in the single- 
period lattice) or both sub-bands (in the doubled-period lattice). On reaching 
the upper-energy state they are partly allowed to leave the lattice towards 
the continuum. For the evaluation of the number of transmitted particles we 
connect the system to incoming and outgoing leads, which mimic its coupling 
to the continuum by injecting and extracting a steady-state particle current. 
A schematic representation of the two-band situation is given in the top part 
of Fig. 2. 




2 3 4 5 



Fig. 2. Top: schematic representation of the two-band system connected to an in- 
coming lead (1) and an outgoing lead (5), showing the half-period Tg/2 of Bloch 
oscillations and the tunnelling time r/2. Bottom: the equivalent birefringent system 
for light propagation. 
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The transmittivity of a BEC of 87 Rb atoms has been evaluated by using the 
scattering matrix formalism adapted to the case of out-of-equilibrium leads [8] . 
We have used the set of parameters U = 3.5 E r with E r = h 2 k 2 /2m, f3 2 ~ 
0.01 and k = 8.2 /mi -1 . In the single-period case the particle current varies 
monotonically with the external force and hence with the period of Bloch 
oscillations. After period doubling the boson wavepacket is split at the edge of 
the Brillouin zone (point (a) in Fig. 2), where it can either be Bragg-reflected 
to point (b) or tunnel into the upper sub-band. The interference between 
the wavepackets reaching point (b) by these various paths gives the outgoing 
current shown in Fig. 3. The minima in transmittivity towards the continuum 
are located at integer values of T b /t, where r = (3ir 2 /8)(hN/n s AE). 
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Fig. 3. Interference pattern in condensate transmittivity from period doubling, as a 
function of Tb/t. 



3.1 Optical equivalent 

A condensate wavepacket propagating in an infinite doubled-period lattice 
is equivalent to a light beam of frequency u travelling through a five-layer 
optical medium (see Fig. 2). The first and the last medium in the bottom part 
of Fig. 2 are semi-infinite and play the role of the two leads. The second and 
the fourth layer stand for the two energy bands, while the middle layer mimics 
the minigap. 

Let ttj and be the transmission and reflection coefficient at the interface 
between the media % and j, connected by tij = 1 F° r a suitable coupling 

between the leads and the lattice t\^ an d £4,5 have unitary modulus and on 
crossing these interfaces the wave takes up an irrelevant phase factor. At the 2- 
3 and 3-4 interfaces we can mimic the effect of the minigap in Bragg-reflecting 
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the BEC by imposing total reflection of light via r 2 ^ = exp(io2,3) and = 
exp(ia 3 4), and allow for tunnel by propagation through layer 3 only via an 
evanescent wave. If the media 2 and 4 have the same refractive index and the 
same optical depth, we have 0:3,4 = a 2 ,3 — n and by symmetry we can set 
0:2,3 = — «3,4 = tt/2. We can then use the recursive formula [15] 



Tid+2 ~ l + r iJ+1 r j+ld+2 e*°^ (6) 

to calculate the transmission coefficient between the media 1 and 5, where Oj 
is the phase shift acquired by light travelling through the j-th layer. To pursue 
the analogy with the doubled-period lattice we have set o 2 = o 4 = uT B /2 and 
03 = iur/2, with r = lixjuj. The total transmission coefficient \t\$\ 2 for the 
light intensity is then found to have minima when the ratio Tb/t takes integer 
values, as in the case of the doubled-period lattice. 



The main difference between the two patterns is that in the BEC case the 
height of the peaks is largest at low values of Tb/t (see Fig. 3), whereas in 
the optical analog all peaks have the same height. Decreasing Tg leads to a 
decrease of the optical depth of media 2 and 4 in the five-layer system, while 
for the condensate it means that the bosons leave the lattice towards the 
continuum after having travelled through a lower number of sites. Therefore, 
decreasing Tg is equivalent in this case to shortening the array and hence helps 
the tunnelling. 



4 Interference from a Fibonacci chain 



The symmetry of an optical potential that is created by the interference of 
optical laser beams is completely determined by the geometric arrangement 
of the beams. Therefore, one can not only realize lattices with various sym- 
metries in the laboratory, but also design quasiperiodic optical potentials. 
Here we give as an example the use of the projection method that takes ori- 
gin from solid state physics for generating a Fibonacci chain (see left panel 
in Fig. 4). The ID Fibonacci arrangement is obtained from a 2D periodic 
square lattice by projecting all sites belonging to a strip with the irrational 
slope o = arctan(2/ (a/5 + 1)) onto a line with the same slope (see for instance 
[16]). With this construction the distance between two neighbouring projected 
sites can take two different values, A or B say. Their ratio A/B = (y / 5 + l)/2 is 
the so-called golden ratio and the sequence of distances follows the Fibonacci 
chain rule ABB AB ABB • • • . This sequence can be obtained by the transfor- 
mation rule A — > B and B — > BA. 



7 



Fig. 4. Left panel: the projection method for creating the Fibonacci chain. Right 
panel: the equivalent set-up for a Fibonacci optical chain realized within five laser 
beams. 



This method can be applied to an atomic gas by using a five-laser configuration 
as is shown in the right panel of Fig. 4. Four laser beams build the 2D square 
lattice by their interference. The fifth laser is aligned with the longitudinal 
axis of the cigar-shaped magnetic trap and drives the condensed bosons along 
the direction of the Fibonacci array. Within this model the hopping energies 
follow the Fibonacci sequence, but in order to make a direct comparison with 
the results on the periodic ID lattices we consider below the case in which the 
site energies (rather than the bond lengths) form the Fibonacci chain. This 
scheme should at least qualitatively give the correct physical picture of matter 
waves propagating through a quasi-periodic array. 



The total density of states for a very long Fibonacci array is shown in the left 
panel of Fig. 5. The fragmentation of the spectrum is typical of quasi-periodic 
and aperiodic systems, as it is well known from previous solid state studies (see 
for instance Ref. [17]). In particular, in the classical case of a quasi-periodic 
Fibonacci chain the spectrum is known to be a Cantor set with measure zero. 
For a chain of 100 wells, the site- projected density of states (see right panel in 
Fig. 5) has a general envelope resembling that of the single-period lattice, but 
is modified by the quasi-periodicity favouring the population of certain sites. 
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Fig. 5. Left panel: total DOS of a very long Fibonacci chain (1000 sites) as a function 
of (E — Eo)/t. Right panel: projected DOS for a BEC driven by a constant force 
through a Fibonacci chain of 100 sites at constant energy E = Eq — 2t, as a function 
of the site number. 



The condensed bosons travelling through the chain explore an energy spec- 
trum which on average is rather more akin to the single-period band structure 
than to the doubled-period one. Nevertheless, the transmittivity through a 
Fibonacci chain presents an interference pattern (see Fig. 6). The presence of 



peaks is a signature of quasi-periodicity leading to 
energy spectrum. 



'quasi-minigaps" in the 
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Fig. 6. Interference pattern in condensate transmittivity through a Fibonacci chain, 
as a function of Tb/t. 
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5 Conclusions 



In summary, we have reviewed earlier work on condensate transport in optical 
lattices and proposed an optical equivalent for the interference pattern that is 
generated by the opening of a minigap in the energy spectrum when matter 
waves are propagating in a lattice with doubled period. 

We have then shown that condensate interference also results from the open- 
ing of sharp depressions in the spectral density of states for a matter wave 
propagating in a quasi-periodic array. Although quasi-periodicity is often said 
to be in some sense intermediate between perfect periodicity and complete 
disorder, it has been shown in earlier work [8] that in neither of these two 
cases an interference pattern of any sort is found in the absence of minigaps. 
We have also proposed a method by which a quasi-periodic modulation of 
the Fibonacci type may be created by optical means for a condensate in an 
elongated magnetic trap. 

Acknowledgements 

This work was partially supported by INFM through the PRA-Photonmatter. 
References 

[1] A. Minguzzi, S. Succi, F. Toschi, M. P. Tosi, P. Vignolo, Phys. Rep. In press. 

[2] M. Kozuma, L. Deng, E. W. Hagley, J. Wen, R. Lutwak, K. Helmerson, S. L. 
Rolston, W. D. Phillips, Phys. Rev. Lett. 82 (1999) 871. 

[3] S. Inouye, T. Pfau, S. Gupta, A. P. Chikkatur, A. Gorlitz, D. E. Pritchard, 
W. Ketterle, Nature 402 (1999) 641. 

[4] M. B. Dahan, E. Peik, J. Reichel, Y. Castin, C. Salomon, Phys. Rev. Lett. 76 
(1996) 4508. 

[5] B. P. Anderson, M. A. Kasevich, Science 282 (1998) 1686. 

[6] F. S. Cataliotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trombettoni, 
A. Smerzi, M. Inguscio, Science 293 (2001) 843. 

[7] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, I. Bloch, Nature 415 (2002) 
39. 

[8] P. Vignolo, Z. Akdeniz, M. P. Tosi, J. Phys. B 36 (2003) 4535. 



10 



[9] J. C. Slater, Phys. Rev. 87 (1952) 807. 

[10] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86 (2001) 2353. 

[11] R. Farchioni, G. Grosso, G. Pastori Parravicini, Phys. Rev. B 53 (1996) 4294. 

[12] C. Fort, F. S. Cataliotti, L. Fallani, F. Ferlaino, P. Maddaloni, M. Inguscio, 
Phys. Rev. Lett. 90 (2003) 140405. 

[13] P. Vignolo, R. Farchioni, G. Grosso, Phys. Rev. B 59 (1999) 16065. 

[14] R. Farchioni, G. Grosso, P. Vignolo, Phys. Rev. B 62 (2000) 12565. 

[15] M. Born, E. Wolf, Principles of Optics, Pergamon, London, 1959. 

[16] T. Fujiwara and T. Ogawa, Quasicrystals, Springer, Berlin, 1990. 

[17] J. Vidal, D. Mouhanna and T. Giamarchi, Phys. Rev. B 65 (2001) 014201. 



11 



